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C^i Abstract 

a.: 

We calculate the independent helicity amplitudes in the decays B — > K*£ + £ and 
-B — ► in the so-called Large-Energy-Effective-Theory (LEET). Taking into ac- 
count the dominant 0(a s ) and SU(3) symmetry-breaking effects, we calculate various 
^ ■ Dalitz distributions in these decays making use of the presently available data and de- 

cay form factors calculated in the QCD sum rule approach. Differential decay rates in 
the dilepton invariant mass and the Forward-Backward asymmetry in B —* K*£ + £~ 
are worked out. We also present the decay amplitudes in the transversity basis which 
has been used in the analysis of data on the resonant decay B — ► K*J/ip(-^> £ + £~). 
Measurements of the ratios Ri{s) = dT Ht (s)(B -> K*£+£ )/dT Hi (s)(B -> p£v e ), in- 
volving the helicity amplitudes i = 0, +1, —1, as precision tests of the standard 
model in semileptonic rare -B-decays are emphasized. We argue that Ro(s) and R-(s) 
can be used to determine the CKM ratio Vt s \ and search for new physics, where 
the later is illustrated by super symmetry. 
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1 Introduction 



Rare B decays involving flavour-changing-neutral-current (FCNC) transitions, such as b — > 
57 and b — > s£ + £~, have received a lot of theoretical interest [1]. First measurements of 
the decay B — > X s 7 were reported by the CLEO collaboration [2]. These decays are now 
being investigated more precisely in experiments at the B factories. The current world 
average based on the improved measurements by the CLEO [3] , ALEPH [4] and BELLE [5] 
collaborations, B[B — > X s ^) = (3.22±0.40) x 1CT 4 , is in good agreement with the estimates 
of the standard model (SM) [6-8] , which we shall take as B(B -> X sl ) = (3.50±0.50) x 10" 4 , 
reflecting the parametric uncertainties dominated by the scheme-dependence of the quark 
masses. The decay B — > X s 7 also provides useful constraints on the parameters of the 
super symmetric theories, which in the context of the minimal super symmetric standard 
model (MSSM) have been recently updated in [9]. 

Exclusive decays involving the b — > s^y transition are best exemplified by the decay B — > 
fT*7, which have been measured with a typical accuracy of ±10%, the current branching 
ratios being [3,10,11] B{B ± -> K*±i) = (3.82 ±0.47) x 10" 5 and B(B° -> K*°~f) = (4.44 ± 
0.35) x 10~ 5 . These decays have been analyzed recently [12-14], by taking into account 
0(a s ) corrections, henceforth referred to as the next-to-leading-order (NLO) estimates, in 
the large-energy-effective-theory (LEET) limit [15,16]. As this framework does not predict 
the decay form factors, which have to be supplied from outside, consistency of NLO-LEET 
estimates with current data constrains the magnetic moment form factor in B — > K*^ in 
the range Tj (0) = 0.27 ± 0.04. These values are somewhat lower than the corresponding 
estimates in the lattice-QCD framework, yielding [17] T x x *(0) = 0.32loo2> anc ^ m ^ ne light 
cone QCD sum rule approach, which give typically Tf* (0) = 0.38 ± 0.05 [18,19]. (Earlier 
lattice-QCD results on B — > K*^ form factors are reviewed in [20].) It is imperative 
to check the consistency of the NLO-LEET estimates, as this would provide a crucial 
test of the ideas on QCD-factorization, which have been formulated in the context of 
non-leptonic exclusive 5-decays [21], but which have also been invoked in the study of 
exclusive radiative and semileptonic 5-decays [12-14]. The decays B — > pj and B — > K*^y 
provide a good consistency check of this framework, with the branching ratios, the isospin- 
violating ratio ^pp°^°7) — ^ an< ^ direct CP-violating asymmetries, such as A(p ±, y) = 
b\b-^p-^]+b\b+^p+'1) ' being the quantities of interest [12,14]. Likewise, isospin-violation in 
the decays B — > K*j, defined as A _ = p|go^^ t o^~p^_^^»_^ and its charge conjugate 
A +, will also test this framework [22]. 

The exclusive decays B — > K*£ + £~, £ ± = e ± ,fj, ± have also been studied in the NLO- 
LEET approach in [13,23]. In this case, the LEET symmetry brings an enormous simplicity, 
reducing the number of independent form factors from seven to only two, corresponding to 
the transverse and longitudinal polarization of the virtual photon in the underlying process 
B — > fT*7*, called hereafter ^*(q 2 ) and ^*{q 2 )- The same symmetry reduces the number 
of independent form factors in the decays B — >• plvi from four to two. Moreover, in the 
g 2 -range where the large energy limit holds, the two set of form factors are equal to each 
other, up to S'[/(3)-breaking corrections, which are already calculated in specific theoretical 
frameworks. Thus, knowing V u b precisely, one can make theoretically robust predictions 
for the rare 5-decay B — > K*£ + £~ from the measured B — > p£vi decay in the SM. The 
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LEET symmetries are broken by QCD interactions and the leading 0(a s ) corrections in 
perturbation theory are known [13,23]. We make use of these theoretical developments and 
go a step further in that we calculate the various independent helicity amplitudes in the 
decays B — > K*£ + £~ and B — > plve in the NLO accuracy in the large energy limit. We 
recall that a decomposition of the final state B — > K*(— > Ktt)£ + £~ in terms of the helicity 
amplitudes H± R (q 2 ) and HQ ,R (q 2 ), without the explicit 0(a s ) corrections, was undertaken 
in a number of papers [24-29]. In particular, Kim et al. [26,27] emphasized the role of 
the azimuthal angle distribution as a precision test of the SM. Following closely the earlier 
analyses, we now calculate the 0(a s ) corrections in the LEET framework. 

Concentrating on the decay B — > K*£ + £~, the main theoretical tool is the factorization 
Ansatz which enables one to relate the form factors in full QCD (called in the literature 
A (q 2 ), A^q 2 ), A 2 (q 2 ),V(q 2 ),T 1 (q 2 ),T 2 (q 2 ),T 3 (q 2 )) and the two LEET form factors £ ± (q 2 ) 
and £||(g 2 ) [13,23]; 

f k (q 2 ) = c^ ± ( q 2 ) + q,£ii(g 2 ) + $ b ® T k ® $ y , (1) 

where the quantities Cj (i =_L, ||) encode the perturbative improvements of the factorized 
part 

Q = Cf ) + ^C\ l) + (2) 

7T 

and Tk is the hard spectator kernel (regulated so as to be free of the end-point singularities), 
representing the non-factorizable perturbative corrections, with the direct product under- 
stood as a convolution of T k with the light-cone distribution amplitudes of the B meson 
(<&b) and the vector meson ($y). With this Ansatz, it is a straightforward exercise to im- 
plement the 0(a s )-improvements in the various helicity amplitudes. The non-perturbative 
information is encoded in the LEET-form factors, which are a priori unknown, and the 
various parameters which enter in the description of the non-factorizing hard spectator 
contribution, which we shall discuss at some length. The normalization of the LEET form 
factor £;f*(g 2 ) at q 2 = is determined by the B — > K*^ decay rate; the other form factor 
^ K (q 2 ) has to be modeled entirely for which we use the light cone QCD sum rules. This in- 
put, which for sure is model-dependent, is being used to illustrate the various distributions 
and should be replaced as more precise data on the decay B — > p£u e becomes available, 
which then can be used directly to determine the form factors £f* (q 2 ) and £^*(<? 2 ), taking 
into account the S'L r (3)-breaking effects. 

Using the effective Hamiltonian approach, and incorporating the perturbative improve- 
ments, we calculate a number of Dalitz distributions, the dilepton invariant mass distri- 
bution for the individual helicity amplitudes (and the sum), and the forward-backward 
asymmetry in B — > K*£ + £~ . As the range of validity of the LEET-based estimates in this 
decay is restricted to the large- E K * region, we shall restrict ourselves to the low s-region 
in the dilepton invariant mass, which for the sake of definiteness is taken as s < 8 GeV 2 . 
We shall also neglect the contributions from the long-distance effects to the final state 
B — > K*£ + £~, arising from the process B — > K*(p,u,4>) — > K*£ + £~, as they are expected 
to be tiny due to the CKM-suppression and the small leptonic branching ratios of the 
vector mesons p,u,cj). To project out the various helicity components experimentally, one 
can use the Dalitz distribution in the dilepton invariant mass (s = q 2 ) and cos#x, where 
9k is the polar angle of the K meson in the rest system of the K* meson measured with 
respect to the helicity axis, i.e., the outgoing direction of the K*. The angular distribution 
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allows to separate the 0-helicity component \H (s)\ 2 = \Hq(s)\ 2 + \H^(s)\ 2 and the sum 
\H + (s)\ 2 + \H„(s)\ 2 . In the SM > and otner beyond-the-SM scenarios considered here which 
have the same operator basis, the component H + (s) is negligibly small. This holds for 
both the left-handed and right-handed projections, H+{s) and H R (s). We show this here 
in the case of the SM. Hence, for all practical considerations, these components can be 
ignored and we concentrate on the H_(s) and Hq(s) components. We show the systematic 
improvements in 0(a s ) and 1/M in H_(s) and H (s) in these decays. Their measurements, 
in conjunction with the decay distributions in B — > p£v^ will serve as precision tests of the 
flavour sector in the SM, yielding |14&|/|14s|, and in searching for possible deviations from 
the SM, exemplified here by super symmetry. 

We also work out the decay amplitudes for B — > K*(—> Kn)£ + £~ in the transversity 
basis [30-32], which has been used by several experimental groups to measure the corre- 
sponding amplitudes for the decay B — > K*(— * K7r)J/ijj(— > £ + £~) [33-36]. These involve 
the complex amplitudes Ao(s), A\\(s) and A±(s). The amplitudes in the transversity and 
helicity bases are simply related [37] and, having worked them out in the helicity basis, it is 
a straightforward numerical exercise to work out the moduli and arguments of the ampli- 
tudes in the transversity basis. Restricting ourselves to low-s region (s < 8 GeV 2 ), we show 
the results using the LEET approach both in the LO and NLO. For illustrative purpose, we 
show the amplitudes in the entire kinematically allowed region in the LO. The LEET-based 
transversity amplitudes for the decay B — > K*£ + £~ are found to be in reasonable agreement 
with their measured counterparts in the resonant decay B — > K*J/ip(—> £ + £~). Measure- 
ment of the short-distance component of these amplitudes coming from B — > K*£ + £~ away 
from s = mj^, in particular in the region Am 2 < s < m^, will test the underlying 
LEET-based framework. 

This paper is organized as follows: In section 2, we define the effective Hamiltonian 
and the matrix element for the decay b — > s£ + £~. In section 3, we discuss the form 
factors in the LEET approach for the decay B — > K*£ + £~, borrowing heavily from the 
literature [13,16,23], give parametrizations for the two remaining form factors £±*(s) and 
£^ (s) and specify other input parameters in our analysis. In section 4, we introduce the 

helicity amplitudes H± R (s) and Hq ,r (s), give the 0(a s )-improved expressions for these 
amplitudes and write down the Dalitz distributions in the set of variables (0, s), (cos^x, s), 
and (cos# + ,s). The quantities \H± R (s)\ 2 are shown as functions of s. Likewise, Dalitz 
distributions in (cos#+, s) are shown for the two dominant components, H (s) and H-(s), 
and adding all three components. We also show the dilepton invariant mass distributions 
for the individual helicity amplitudes, and their sum, and the forward backward asymmetry, 
making explicit the 0(a s ) improvements. Section 5 describes the amplitude decomposition 
for B — > K*£ + £~ in the transversity basis. We show the amplitudes | ^4_o (^) 1 2 ; l^||( s )| 2 
and |^4±(s)| 2 , as well as the relative phases <f>\\(s) and (j>±(s), making explicit the 0(a s ) 
improvements in these quantities. Extrapolating the LO results for these quantities to 
the J/ip mass, we compare them with data on B — > K*J/ip(—^ £ + £~). In section 6, we 
turn to the decay distributions in the decay B — > piv^ and display the various helicity 
components, Dalitz distributions, and the dilepton {vd) invariant mass. Estimates of 
the B — > p LEET form factors £^(s) and £m(s), which are scaled from their B — > K* 
counterparts incorporating SU(3)-breaking, are also displayed here. Section 7 is devoted 
to the determination of the ratio of the CKM matrix elements |V^&|/|Vf S | from the ratio 
of the dilepton mass spectra in B — > p£vi and B — > K*£ + £~ decays involving definite 
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dV^K* l + l /ds 

helicity states. In particular, we show the dependence of the ratio R-{s) = dT B^ P wj^~ s — 

and Rq(s), involving the helicity-0 components, on the CKM matrix elements |K&|/|V*s|- 
Section 8 is devoted to an analysis of the ratios Rq(s) and R-(s) to probe for new physics in 
the decay B — > K*£ + £~, and illustrate this using some specific supersymmetric scenarios. 
Finally, section 9 contains a summary and some concluding remarks. 



2 Effective Hamiltonian for b — » s£ + £ 

At the quark level, the rare semileptonic decay b — > s£ + £~ can be described in terms 
of the effective Hamiltonian obtained by integrating out the top quark and W ± bosons: 

Heff = » ( 3 ) 

* i=l 

where Vij are the CKM matrix elements [38] and Gf is the Fermi coupling constant. We 
use the operator basis introduced in [6] for the operators Oi, i — 1, ... ,6, and define: 



o 7 = 



g em m b 
8tt 2 



5^(1+75) bF^, 



g s m b 

' 8tt 2 



0. 



9,10 



2n 



(££)v,a(s b) v _ A , 



s^il + ^mbjG* (4) 



(5) 



where a em = g 2 m /An is the electromagnetic fine-structure constant. T a , a = 1....8 are the 
generators of QCD, and i, j are SU (3) color indices. Here F^ and G^ u denote the electro- 
magnetic and chromomagnetic field strength tensor, respectively. The above Hamiltonian 
leads to the following free quark decay amplitude: 



M(b -> s£ + T) 



GpCien 

V2n 



V*V tb [C 9 [s^Lb] [l-fi] + C w [s^Lb] [l^d] 



-2m b C° 



eff 



si(Juu — Rb 



\i-ft] 



(6) 



Here, L/R = (1 =F7 5 )/2, s = q 2 , = \ [7^,7,,] and <^ = (p+ where p± are the 

four-momenta of the leptons. We put m s /m b = and the hat denotes normalization in 
terms of the 5-meson mass, mg, e.g. s = s/m 2 B , rh b = m b /mB- Here and in the remainder 
of this work we shall denote by m b = m b (n) the MS mass evaluated at a scale //, and by 
m^poie the pole mass of the 6-quark. To next-to-leading order the pole and MS masses are 
related by 



m b (fx) = m b>po i e 



An 



3 In 



+ 0(a 2 s ) 



(7) 



Since we are including the next-to-leading corrections into our analysis, we will take the 
Wilson coefficients in next-to-leading-logarithmic order (NLL) given in Table I. 
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Ci 


c 2 


c 3 


C 4 


<?5 


Ce 


LL 
NLL 


-0.257 
-0.151 


1.112 
1.059 


0.012 
0.012 


-0.026 
-0.034 


0.008 
0.010 


-0.033 
-0.040 




Cf 




c 9 


Cio 




°10 


LL 
NLL 


-0.314 
-0.308 


-0.149 
-0.169 


2.007 
4.154 




-4.261 


4.214 


-4.312 



Table 1: Wilson coefficients at the scale // = 4.6 GeV in leading-logarithmic 
(LL) and next-to-leading-logarithmic order (NLL) [13]. 

3 Form factors in the Large Energy Effective Theory 

Exclusive decays B — > K*£ + l~ are described by the matrix elements of the quark 
operators in Eq. (6) over meson states, which can be parameterized in terms of form 
factors. 

For the vector meson K* with polarization vector e M , the semileptonic form factors of 
the V — A current are defined as 

(K*(p,e*)\(V - A) lx \B{p B )) = -i e;(m B + m K ,)A 1 (s)+t (p B + p)^e*p B ) Us) 



m B + vrtK* 



q^PB) — (M*) - Ms)) + W ^V B f 2V } S) ■ (8) 
Note the exact relations: 



a i \ m B + m K * m B - m K » 

A o (0) = A 3 (0), 

(K*\d li A»\B) = 2m K *(e*p B )A (s). (9) 

The second relation in (9) ensures that there is no kinematical singularity in the matrix 
element at s = 0. The decay B — > K*i + i~ is described by the above semileptonic form 
factors and the following penguin form factors: 

(K*( P ,e*)\Cfrsa^(l + l5 )b\B(p B )) = i e^ p<T e» CT 2T X ( S ) (10) 

+ T 2(s){e*(m| - m 2 K ,) - (e*p B ) (p B + P )^ 

+ T 3 (s)(e*p B ) \q»~ 2 S 2 (PB +p)A ■ 

I rn B ,n K* ) 

The matrix element decomposition is defined such that the leading order contribution from 
the electromagnetic dipole operator O-j reads %{s) = Cj S Ti(s) + . . ., where Tj(s) denote 
the tensor form factors. Including also the four-quark operators (but neglecting for the 
moment annihilation contributions), the leading logarithmic expressions are [43] 



(11) 



T 2 (s) = CfT 2 {s) + Y{s) S A^s), (12) 

2mb{i7iB — rriK*) 



%(s) = C 7 ett T 3 (s) + Y(s) 



m B - m K * m B + m K * 

2m b MS) - 2m h Ms) 



(13) 



with Cf = C 7 - C 3 /3 - 4C 4 /9 - 20^5/3 - 80C 6 /9 = C 7 - (4C 3 - C 5 )/9 - (4C 4 - C 6 )/3, 
and 

Y(s) = h(s,m c )(3C 1 + C 2 + 3C 3 + C 4 + 3C 5 + C 6 ) 

- X - h(s, m b ) (4 (C 3 + C A ) + 3^5 + C 6 ) - X - h(s, 0) (C3 + 3C4) 

+ ^ f^3 + 2C 4 + ^ , (14) 



9 V3 3 

where the "barred" coefficients (7j ( for i=l,...,6) are defined as certain linear combinations 
of the Ci, such that the Q coincide leading logarithmic order with the Wilson coefficients 
in the standard basis [44]. Following Ref. [13], they are expressed as : 

C2 = C2 — — C\ , 
6 

1 8 
C3 — C3 — - C 4 + 1G C 5 — - C 6 , 
o 3 

C 4 = -C4 + 8C6, 

C*5 = C3 — - C4 + 4C 5 — - C e , 
o 3 

C'e = ic 4 + 2C 6 . (15) 



The function 



2 



, arctan . z > 1 , 

4 A m o 2 \ 4 . , n y/z-1 



4 / Tfi 2 \ 4 
h(s,m q ) = --\h 1 -±--- z y-(2 + z) v^lT 



1 + 



27T 



ln ^i T 



(16) 

is related to the basic fermion loop. (Here z is defined as Arrii/ s.) Y(s) is given in the NDR 
scheme with anticommuting 75 and with respect to the operator basis of [6]. Since C 9 is 
basis-dependent starting from next-to-leading logarithmic order, the terms not proportional 
to h(s, m q ) differ from those given in [44]. The contributions from the four-quark operators 
Oi_6 are usually combined with the coefficient C 9 into an "effective" (basis- and scheme- 
independent) Wilson coefficient C 9 cff (s) = C 9 + Y(s). 

Recently, it has been shown that the symmetries emerging in the large energy limit [16] 
relate the otherwise independent form factors entering in the decays of B mesons into light 
mesons. However this symmetry is restricted to the kinematic region in which the energy 
of the final state meson scales with the heavy quark mass. For the B — > K*£ + £~ decay, 
this region is identified as s ~ 8 GeV 2 . 
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Thus, in the large energy limit, the standard form factors A , Ai, A 2 , V, 7\, T 2 and T 3 can 
be expressed in terms of two universal functions £||( s ) and £j_(s) [16] : 



A (s 

A 1 (s 
A 2 (s 
V(s 

T 2 (s 
Us) 



rn 



v 



m B E v 
2E V 



m B ■ 
(1 + 



I- my 
my 

m B 



Us) , 



m B 



\us)-^-m 

iii' - 



V 



(1 + ^) Us) 
m B 

Us) , 
(1 



m, 



m 



-) Us) 



V 



is 



( s ) _ my 

Ey 



mi 



m, 



where 



£V = !^(i__V + 4), 

2 V mi mi/ 



(17) 

(18) 
(19) 

(20) 

(21) 
(22) 

(23) 
(24) 



refers to the energy of the final vector meson V and ^^(s) refer to the form factors in 
the large energy limit (called subsequently as the LEET form factors). However, these 
symmetries are broken by factorizable and non-factorizable QCD corrections, worked out 
in the present context by Beneke et al. [13,23]. Since, we are using in our analysis the 
definitions of the form factors £mi(s) by Charles et al. [16], the factorizable corrections 
obtained in [23] are expressed as follows: 



A 1 (s) = 
Ms) = 



2E 



v 



m B + m v 

m B 
m B — m v 



Us) 
Us) 



47T 

my 

Ey 



■AA 1 , 

*||(*)(l + ^[-2 + 2L]) 



H -r^~ AA 2 

47T 



(25) 
(26) 



with 



L = 



2E> 



v 



In 



2Ev 



AAi = 0, AA 2 
AFu 



m B - 2E V m B 
my m B {m B — 2Ey) 



m B — m v 

.2. 



AEy 



AFu 



v 



N c m B 



(27) 



where f B and fy are, respectively, the meson decay constants for the B meson and the 
corresponding V meson. The above expression for AA 2 also involves a non-perturbative 



quantity (£ + 1 )+. Formally, (£ + 1 )- 



-, but nothing more is known about this universal 



-+ /+• v> V+ /+ " A QCD - 

parameter at present. It is estimated to lie in the range (0.2 — 0.5 GeV)^ 1 [13], following 
which we take (0.3 GeV)^ 1 as our default value for this quantity in our calculations. For 
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the K* meson, we use the result quoted in Ref. [13] : < u^ 1 >\\= 3.75. Concerning the 
form factors A (s) and V(s), defined respectively in Eqs. (17) and (20), they hold exactly 
to all orders in perturbations theory and this defines the factorization scheme. 

The remaining contributions arising from the hard spectator corrections for the B — > 
V l + l~ decay have been computed recently by Beneke et al. [13], yielding 

Us) = 71(a), (28) 

T 2 (s) = —^T ± (s), (29) 
m B 

r 3 ( s )-^-r 2 ( a ) = 7||( s ), (so) 



with 



4"7T 



l\ c m B ±J u Jo 



T _ t m v ( r (o) (*sCf r (i) 



+ ^^ S llE/ ^^bMJ/u^^T^M. (31) 

Here Cf = 4/3, A^ c = 3, Sj_ = 1, = my/Ey, and the hard-scattering term T a ,±(w, a;) 
(a =±, ||) is expanded as : 

Ta,±(u,u) = Ti°i(«,a;) + ^tJ, 1 »(u, W ) , (32) 

where denotes the usual X* decay constant /r-* and fx*,± refers to the (scale- 

dependent) transverse decay constant defined by the matrix element of the tensor current. 

The coefficient (a =_L, ||) in (31) represents the next-to-leading order form factor 
correction, and can be expressed as : 

C?> = + Cf f) , (33) 

(f ) 

where C a contains a factorizable term from expressing the full QCD form factors in terms 
of £ a in Eqs. (11), (12) and (13). The non-factorizable correction cf' 1 is obtained by com- 
puting matrix elements of four-quark operators and the chromomagnetic dipole operator. 
The matrix elements of four-quark operators require the calculation of two-loop diagrams, 
and the result for the current-current operators as well as the matrix element of the 
chromomagnetic dipole operator can be extracted from Ref. [45]. The 2-loop matrix el- 
ements of the QCD penguin operators have not yet been computed and hence will be 
neglected. This should be a very good approximation due to the small Wilson Coefficients 
of the penguin operators. For the definitions of the parameters in Eq. (31), we refer to [13]. 

Lacking a complete solution of non-perturbative QCD, one has to rely on certain ap- 
proximate methods to calculate the above form factors. In this paper, we take the ones 
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A 1 A 2 A V 


Ti T 2 T 3 


F(0) 

Cl 

c 2 


0.294 0.246 0.412 0.399 
0.656 1.237 1.543 1.537 
0.456 0.822 0.954 1.123 


0.334 0.334 0.234 
1.575 0.562 1.230 
1.140 0.481 1.089 



Table 2: Input values for the parameterization (34) of the B — > K* form factors. Renor- 
malization scale for the penguin form factors % is /i — m b [18]. 



M w 


80.4 GeV 


!b 


200 ± 20 MeV 


m t {m t ) 


167 ±5 GeV 


/ft - *, ii 


225 ± 30 MeV 


m b)PO i e (2 GeV) 


4.6 ±0.1 GeV 


/A-!±(lGeV) 


185 ± 10 MeV 


m c 


1.4 ±0.2 GeV 


f P (1 GeV) 


198 ± 7 MeV 




1/137 




(3 ± 1) GeV 


T B 


1.65 ps 




0.2 ±0.1 


\v; s v tb \ 


0.041 ±0.003 




0.05 ±0.1 


R b = \v ub \/\v* s v tb \ 


0.094 ±0.014 


er } (o) 


0.28 ±0.04 


A («/=5) 
QCD 


220 ± 40 MeV 


d p) (o) 


0.22 ±0.04 




0.3 ±0.2 {GeV)- 1 




3.48 



Table 3: Input parameters and their uncertainties used in the calculations of 
the decay rates for B — > K*£ + £~ and B — > pivg in the LEET approach. 

given in [18], obtained in the framework of Light-cone QCD sum rules, and parametrized 
as follows: 

F(s) = F(0)exp( Cl s±c 2 s 2 ). (34) 

The coefficients in this parametrization are listed in Table 2, and the corresponding LEET 
form factors £j_(s) and £ii(s) are plotted in Fig. 1. The range £j_(s) = 0.28 ± 0.04 is 
determined by the B — > K*^ decay rate, calculated in the LEET approach in next-to- 
leading order [12-14] and current data. This gives somewhat smaller values for Ti(0) and 
T 2 (0) than the ones estimated with the QCD sum rules. 



4 Distributions in the Decay B — > K*£ + £ 

We introduce the helicity amplitudes for the decay B — > > i\T(pif)±7r(p w ))£ + (p + )£~(p_), 

which can be expressed as [26]: 

= ± c L , R y/X) , 

H (s) = &l,r 7= + 7= , 35 

myvs my^Js 

where P.L = — — s)/2 and A = [f (m^ — m v — s) 2 — m v s] and V stands here 
for the vector meson K* . Our definitions for the quantities a^R, bL,R and cl,r differ from 
those used by Kim et al. [26] by a factor of 1/y/s. They read as follows: 
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s) 




[BFS] 



1 2 3 4 5 6 7 

s (GeV 2 ) 




1 2 3 4 5 6 7 

s (GeV 2 ) 



0.7 
0.6 
0.5 
0.4 
0.3 




1 2 3 4 5 6 7 

s (GeV 2 ) 



1.1 

1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 




1 2 3 4 5 6 7 

s (GeV 2 ) 



Figure 1: LEET form factors (s) for B — * K* / + / . The two columns denoted by 

[AS] and [BFS] represent, respectively, our ^^j, \s) and the ones used by Beneke et ai in 
ref [13]. The central values are represented by the dashed curves, while the bands reflect 
the uncertainties on the form factors. 
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dL/R = 
h L/R = 



i(m B + m v ) 
2 m b m B y[s~ 

i 



m b m B (ra| - m v ) y/s 



s m B (±C 10 - C 9 )^i(s) + 471 (s) m b (m v - m B )E v ^ , 
4 Ti(s) m 6 (-m| + m^) 



(36) 
(37) 



+m B s( - 2 m fe |ri(s) + T 3 (s) - T 2 (s)| + A 2 (s) (±C 10 - C 9 )(m B - m v )) 



271 (s) mi, (m B + m v ) + (tC w + C 9 ) s V(s) 



m b (m B + my) y/s 

We show the helicity amplitudes \H^(s)\ 2 , \H^(s)\ 2 , \H R (s)| 2 , and \H R (s)\ 2 in Fig. 2, Fig. 3, 
Fig. 4, and Fig. 5, respectively. 



(38) 



4.1 Dalitz distributions 

Using the above helicity amplitudes, the angular distribution in B — > > Kn)£ + t 
is given by the following expression: 



d 4 T 



3 a^GVAmglWEl 2 



ds d cos 9 K d cos 9 + d(f) 128(27r) 6 m| 
x{4 cos 2 ^ sin 2 # + (|/7 R (s)| 2 + |/7 L (s)' 2 ' 

+ sm 2 6 K (l + cos 2 #+) {\H^(s)\ 2 + \H^s)\ 2 + \H R (s)\ 2 + \H R (s)\ 2 ^ 

-2 sin 2 9 K sin 2 6 + cos20 Re^(s)/7 R *(s) + H^(s)H L *(s)^ 

-sin 20 Im(/7^(s)/7 R *(s) + J7^(s)/7 L *(s)) 



(39) 



— sin 2^ sin2# H 



cos0 Re{H R (s)H R *(s) + H R (s)H R *(s) + + H L _(s)H^{s) 



- sin0 Im(i?f(s)<( S ) - H R {s)H R \s) + iff - H^H^a] 

-2 sin 2 ^ cos# + (|77f(s)| 2 - |77^(s)| 2 - |77^(s)| 2 + |77f(s)| 2 ) 



+2 sin# + sin 20* cos0 Re^(s)^(s) - H R (s)H R *(s) - H^(s)H^*(s) + H L (s)H^*(s] 
-sin<l>lm(H*(8)Ht(8) + H?( 8 )Ht(s) - H^s)H L *(s) - H L (s)H L *(s))] } . 

Here, the various angles are defined as follows: 9 K is the polar angle of the K meson in the 
rest system of the K* meson, measured with respect to the helicity axis, i.e., the outgoing 
direction of the K*. Similarly, 9 + is the polar angle of the positively charged lepton £ + in 
the dilepton rest system, measured with respect to the helicity axis of the dilepton, and is 
the azimuthal angle between the two planes defined by the momenta of the decay products 
K* -> Kn and 7* -> 

Integrating over the angle 6k and 9 + , we get the Dalitz distribution in the remaining 
two variables (0, s): 

~ 2 ^v^^|WJ 2 ^{|/7 ( S )| 2 + |/7 + ( S )| 2 +|/7_( S )| 2 (40) 



d 2 B a 2 em G 2 F 
3 384tt 5 



d(j) ds 
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-cos 20 Re(H*(s)H*(s) + HL(s)Ht*(s)) + sin 20 Im (h r (s)H r *(s) + H^(s)H^(s)j } 
where r B is the 5-meson life time, and : 



\H 


»| 2 = 




» 


I 2 + \H* 


s)\ 2 , 




\H+ 


»r = 




» 


l 2 +l^ 


:*)i 2 , 




\H- 


»r = 




» 


\ 2 + \H R 


»r • 


(41) 



Similarly, we can get the Dalitz distributions in (9 k, s) and (9 + , s), which read as follows: 

+ sin 2 ^ (\H + (s)\ 2 + \H_(s)\ 2 )} . (42) 



d cos 9k ds 



rf 2 B aLGi 



dcos9+ ds 



+ (l + cos# + ) 2 |^(s)| 2 + (l-cos^ + ) 2 |^(s)| 2 
+ (1 - cos# + ) 2 |#f(s)| 2 + (1 + cos# + ) 2 \H R (s)\ 2 } 
d 2 B\ Hi) \2 d 2 B\ H _\2 d 2 B\ H+ \ 



dcos6> + (is <icos# + <is dcos9 + ds 

In Figs. (6), (7) and (8), we plot, respectively, the Dalitz distribution given by the two 
dominant partial contributions and the complete expression given in Eq. (43). 



4.2 Dilepton invariant mass spectrum 

The dilepton invariant mass spectrum can be obtained by integrating over the angle 
variables, yielding: 

dB a 2 G 2 - 2 



^^f^lW| 2 {l^( s )| 2 +l^-( s )| 2 + l^)r} • (45) 
In LEET, the helicity amplitudes (35) are expressed as: 



H L J R (s) 



4 Ti(s) m b (m B — m v ) (m B + m v ) Ey 



2 m b m B (m B + m v )^fs 
+ (±Cio - C 9 ) m B (m B + m v ) 2 s A t (s) 

+2t71bV / a{2T 1 (s) m h (m B + m v ) + frC w + C 9 ) s V(s)}}} , (46) 



+ (±Ci - C 9 )m B (m B + my) 2 s A ± (s) 

-2m B v / A{2T 1 (s) m 6 (m B + m v ) + (=FCio + C 9 )s\/(s)}}] , (47) 
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8 (GeV 2 ) 

Figure 2: The helicity amplitude |if+(s)| 2 at next-to- leading order (solid center line) and 
leading order (dashed). The band reflects theoretical uncertainties from the input param- 
eters. 
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Figure 3: The helicity amplitude \H^(s)\ 2 at next-to- leading order (solid center line) and 
leading order (dashed). The band reflects theoretical uncertainties from the input param- 
eters. 
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0.4- 




Figure 4: The helicity amplitude \H+(s)\ 2 at next-to-leading order (solid center line) and 
leading order (dashed). The band reflects theoretical uncertainties from the input param- 
eters. 
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250 
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Figure 5: The helicity amplitude \H^(s)\ 2 at next-to-leading order (solid center line) and 
leading order (dashed). The band reflects the theoretical uncertainties from the input 
parameters. 
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H o /R ( s ) = i 7 rr 8 \m b T 1 (s)\2{m B - m 2 v )E v + m B s\ 

A m b m B m v (—m B + m v )s I I J 

+4 A m B s{2 m b (%(s) - ^j- T 2 (s)) - A 2 (s)(±C w - C 9 )(m B - m y )} 

+ (m B - m v )(m B + m v ) 2 (m B - m v - s)|4 T x {s) m b E v (-m B + m v ) 

+sm B Ai(s)(±C7 10 -C 9 )}] . (48) 

In Figs. (9), (10) and (11) we have plotted, respectively, the dilepton invariant mass spec- 
trum dB\H_\ 2 /ds, dB\H \ 2 /ds and the total dilepton invariant mass, showing in each case the 
leading order and the next-to-leading order results. The contribution proportional to the 
helicity amplitude H+(s) is negligible, and hence not shown, but is is included in calculating 
the total dilepton spectrum. As can be seen from Figs. (9) and (11) the total decay rate 
is dominated by the contribution from the helicity component. The next-to-leading 
order correction to the lepton invariant mass spectrum in B — > K*£ + £~ is significant in the 
low dilepton mass region (s < 2 GeV 2 ), but small beyond that shown for the anticipated 
validity of the LEET theory (s < 8 GeV 2 ). Theoretical uncertainty in our prediction is 
mainly due to the form factors, and to a lesser extent due to the parameters A^ + and the 
£>-decay constant, f B . 
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dB\ H _\2/ds 10" 7 




1 2 3 4 5 6 7 

s (GeV 2 ) 



Figure 9: The dilepton invariant mass distribution dB\n_\ 2 /ds for B — > K*£ + £~ at next- 
to- leading order (solid center line) and leading order (dashed). The band reflects the 
theoretical uncertainties from input parameters. 



0.5; 
0.4 : 



dB Wa? /ds 10- 7 0.3; 

0.2 : 




0.1- 



1 2 3 4 5 6 7 

s (GeV 2 ) 

Figure 10: The dilepton invariant mass distribution dB\H \ 2 /ds for B — > K*£ + £~ at next-to- 
leading order (solid center line) and leading order (dashed). The band reflects theoretical 
uncertainties from input parameters. 
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dB/ds 10" 7 




1 2 3 4 5 6 7 

s (GeV 2 ) 



Figure 11: The dilepton invariant mass distribution for B — ► K*£ + £~ at next-to-leading 
order (solid center line) and leading order (dashed). The band reflects theoretical uncer- 
tainties from the input parameters. 



4.3 Forward-Backward asymmetry 

The differential forward-backward asymmetry (FBA) is defined as [46] 

cLAfb = _ r [S) d{L ^_ + f° dil ^L (49) 
ds J Q duds J-u(s) duds 

The kinematic variables (s, u) are defined as follows 



n2 ' 



which are bounded as 



(50) 

m B 

u = (Pb -P-f - (Pb -P+f (51) 



(2m,) 2 < s <{l-m K *)\ (52) 
—u{s) < u < u(s) , (53) 



with rhf = me/mB, and 



m 



B 



^) = ^\/A(l-4^). (54) 



Note that the variable u corresponds to 8 + , the angle between the momentum of the B- 
meson and the positively charged lepton £ + in the dilepton CMS frame through the relation 
u = —u(s) cos6L|_ [46]. 
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Figure 12: Forward-backward asymmetry cIAfb(B —>■ K*l + l~) / ds at next-to-leading order 
(solid center line) and leading order (dashed). The band reflects the theoretical uncertain- 
ties from the input parameters. 



At the leading order, the FBA in B — > K*£ + £ decays reads as follows 



dA 



FB 



G% a? m 



F "em " b B 



ds 



2% 5 
xC w 



V*V tb \ 2 su(s) 



(55) 



-C 7 cS )^(-l + rnl,+s)+2 ^(C 7 cS % + Re[C 9 cS ]) 
s ' rriB s 



The position of the zero of this function, Sq, is given by solving the following equation: 



Re(C 9 ett ( So )) = -^C 7 
so 



™b„ eff {l-rh 2 K* - s 



1 + rh 2 K * — s 



(56) 



Our results for FBA are shown in Fig. 12 in the LO and NLO accuracy. We essentially 
confirm the results obtained in the NLO-LEET context by Beneke et al. [13]. 



5 Transversity Amplitudes for B — > K*£ + £ and Com- 
parison with Data on B — > K*J/ip(^ £ + £~) 

The decay B — > J/ipK* is described by three amplitudes (Ai]i = 0, ||,_L) in the 
transversity basis, where Ao(s), A\\(s) and ^4j_(s) have CP eigenvalues +1, +1 and — 1, 
respectively [30,32], and should not be confused with the form factors A (s), Ai(s) etc. 
Here, Ao(s) corresponds to the longitudinal polarization of the vector meson K* and A\\(s) 
and A±(s) correspond to parallel and transverse polarizations, respectively. The relative 
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phase between the parallel (transverse) amplitude and the longitudinal amplitude is given 

by (f)\\(s) = aig(^A\\(s) / Ao(s)^ (^(f)j_(s) = &rg(^A±(s) / Ao(s)J j . The transversity frame 

is defined as the J/ip rest frame (see Fig. 13). The K* direction defines the negative x 
axis. The Kit decay plane defines the (x,y) plane, with y oriented such that p y (K) > 0. 
The z axis is the normal to this plane, and the coordinate system is right-handed. The 
transversity angles 9 tr and (p tr are defined as the polar and azimuthal angles of the positively 
charged lepton from the J/vfj decay; 9k* is the K* helicity angle defined in the K* rest 
frame as the angle between the K direction and the direction opposite to the J/ip. This 
basis has been used by the CLEO [33], CDF [34], BABAR [35], and the BELLE [36] 
collaborations to project out the amplitudes in the decay B — > J/ipK* with well-defined 
CP eigenvalues in their measurements of the quantity sin 2(3, where f3 is an inner angle of 
the unitarity triangle. We also adopt this basis and analyze the various amplitudes from 
the non-resonant (equivalently short-distance) decay B K*£ + £~ . In this basis, both the 
resonant B — > K*J/ip — > K*£ + £~ (already measured) and the non-resonant (B — > K*£ + £~) 
amplitudes turn out to be very similar, as we show here. 

The angular distribution is given in terms of the linear polarization basis (A±i(s) = 
(A {l (s) ± A±(s))/V2) and A(s) by 



d 4 T 



ds dcos9 tr d cos 9 K * d4> tr 



= AH • \A (s)\ 2 + f 2 (w) ■ lA^s)] 2 + f 3 (w) ■ \A ± (s)f 

+ r)f 4 (w) ■ lm(Al(s)A ± (s)) + f 5 (w) ■ ReW)A\\(s)) 
+ r)f 6 (w)-lm(A* (s)A±{s)) , 



where r] = +1(— 1) for B° and B + (B° and B ), and the coefficients /i=i,...,6, which depend 
on the transversity angles w = (9k*, 9tr, <f>tr), are given by: 



fi(w 

f'2{w 

h{w 
h{w 
fe(w 



9/(32tt) • 2cos 2 ^ x *(l - sin 2 ^ r cos 2 tr ), 
9/(32tt) • sin 2 9 K *(l - sin 2 9 tr sin 2 tr ), 
9/(32tt) • sin 2 9 K * sin 2 9 tr , 
9/ (327r) • sin 2 9k* sin 29 tr sin <f) tr , 
9/(32tt) • l/v / 2-sin2^sin 2 ^ r sin20 tr , 
9/(32tt) • 1/V2 ■ sin26»^. sin 26^ cos tr . 



In terms of the helicity amplitudes H±[ , introduced earlier, the amplitudes in the linear 
polarization basis, Aq j_ 11 , can be calculated from the relation: 



A +1 (s) = K (Hi(s) + H*(sj) , 
A-!(s) = K(Ht(s) + H?(sj) , 



with k 2 



384vr 5 



Va 



m%\ v tb v ts\ 



Experimental results are conventionally expressed in terms of the spin amplitudes ^4o,±,|| 
normalized to unity, with \A \ 2 + |^4±| 2 + |^4|i| 2 = 1- We show the polarization fractions, 
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Figure 13: Definitions of the transversity angles 9 tr , <\>tn and 9k*- The angles 9 tr and (j) tr are 
determined in the J/ip rest frame. The angle 9k* is determined in the K* rest frame. 

r /r = |.4q(s)| 2 , r*]|/r = |*4.||(s)| 2 and r_|_/T = |-A_i_(.s)| 2 in the leading and next-to-leading 
order for the decay B — > K*£ + £~ in Figs (14), (15) and (16), respectively. Since the 
interference terms in the angular distribution are limited to Re(.4||.4.Q), 1^h{-A±Aq) and 
Im(^4^^4*|), there exists a phase ambiguity: 

0|| -» -0|| , (57) 
4>± -> ±?r - (j>i_ , (58) 
± - 0n _> ±7r - (0 ± - 0n) . (59) 

To avoid this, we have plotted in Figs. (17) and (18) the functions cos0m(s), sm^n(s), 
and cos0^(s), sin0j_(s), respectively, showing their behaviour at the leading and next-to- 
leading order. The dashed lines in these figures correspond to using the LO amplitudes, 
calculated in the LEET approach. In this order, the bulk of the parametric uncertainty 
resulting from the form factors cancels. Although, strictly speaking, the domain of validity 
of the LEET-based distributions is limited by the requirement of large energy of the K* 
(which we have translated into approximately s < 8 GeV 2 ), we show this distribution for 
the entire s-region allowed kinematically in B — > K*£ + £~. The shaded curves correspond to 
using the NLO contributions in the LEET approach. We compare the resulting amplitudes 
|-4.o| 2 ) l-4±| 2 ; l^-||| 2 5 0||( s )> an d <f>±{s) at the value s = m 2 j/^ with the corresponding results 
from the four experiments in Table 4. In comparing these results for the phases, we had to 
make a choice between the two phase conventions shown in Eq. (59) and the phases shown 
in the last row of this table correspond to adopting the lower signs in these equations. We 
note that the short-distance amplitudes from the decay B —* K*£ + £~ are similar to their 
resonant counterparts measured in the decay B — > J/ipK*. We also note that a helicity 
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analysis of the decay B — > J/vpK* has been performed in the QCD factorization approach 
by Cheng et al. [47]. 

The structures in the phases shown in Fig. (17) and Fig. (18) deserve a closer look. 
We note that at the leading order, the phases 0j_(s) and (f>\\(s) are given by the following 
expressions: 



Ms) 



Arg 



Arg 



m b m B 
-tE v U(s) 



-Arg[A(s)], 



-^{sC 9 cS + 2m b m B C 7 cS }Us)} - Arg[A (s)} , 
s C 9 cS + 2m b m B C 7 cff ) - 2 m b m B (c 7 eS + 



2 m b m B 



-Y(s 



,m 



where we can neglect the term proportional to (my/m B ) in the latter equation. The phase 
o (s) = Arg[„4o( s )] is constant in the entire phase space, as shown in Fig. (19). The 
functions in the square brackets in Eqs. (60) and (61) are purely imaginary. However, 
due to the fact that in the SM the coefficients Cg ff and Cf 1 have opposite signs, these 
phases become zero at a definite value of s, beyond which they change sign, yielding a step- 
function behaviour, shown by the dotted curves in the functions cos0||(s) and cos0j_(s) in 
Fig. (17) and Fig. (18), respectively. The position of the zero of the two functions, denoted, 
respectively, by Sq and sjj, are given by solving the following equations: 



Arg 



_L 



m b m B v Sq 
i E v 



{4 C 9 eff (4) + 2 m b m B C 7 cff } U4) 
Arg pl^ {4c^) + 2m b m B C 7 cS } ^ (4 



m b y s ' 



M4) , 

00 (Sq). 



(62) 
(63) 



For the assumed values of the Wislon coefficients and other parameters, the zeroes of the 
two functions, namely s[j and Sq, occur at around s ~ 3 GeV 2 , in the lowest order, as 
can be seen in Figs. (17) and (18), respectively. The LO contributions in sin 0j| (s) and 
sin0j_(s) are constant, with a value around 0, with a small structure around s ~ 3 GeV 2 , 
reflecting the sign flip of the imaginary part in «4m(s) (A±(s)). At the NLO, the phases are 
influenced by the explicit 0(a s ) contributions from the factorizable and non-factorizable 
QCD corrections (see section 3), which also bring in parametric uncertainties with them. 
The most important effect is that the zeroes of the phases as shown for cos0j_(s) and 
cos0||(s) are shifted to the right, and the step-function type bahaviour of these phases in 
the LO gets a non-trivial shape. Note that in both figures a shoulder around s ~ 8 GeV 2 
reflects charm production whose threshold lies at s = Am 2 . 
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Group 


IAI 2 


\M 2 


l^lll 2 


<t>± 


<f>\\ 


CLEO [33] 
CDF [34] 

BaBar [35] 
Belle [36] 

This Work 


0.52 ±0.08 
0.59 ±0.06 
0.60 ±0.04 
0.60 ±0.05 
0.51 


0.16 ±0.09 

o n+ - 13 

u - io -0.11 

0.16 ±0.03 
0.19 ±0.06 
0.21 


0.32 ±0.12 
0.28 ±0.12 
0.24 ±0.04 
0.21 ±0.08 
0.28 


-3.03 ±0.46 
-2.58 ±0.54 
-2.97 ±0.17 
-3.15 ±0.21 
-3.25 


-3.00 ±0.37 
-2.20 ±0.47 
-2.50 ±0.22 
-2.86 ±0.25 
-3.04 



Table 4: Current measurements of the decay amplitudes in the transversity basis for the 
decay B — > J/ipK* . The corresponding amplitudes for the non-resonant decay B — > 
K*l + £~ worked out in this paper in the LO approximation at m 2 t+t _ = ?n 2 j/^ are given in 
the last row. 
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Figure 14: The helicity amplitude |^lo(<s)| 2 in 5 — > K*£ + £~ at next-to-leading order (center 
line) and leading order (dashed). The band for NLO reflects theoretical uncertainties from 
input parameters. 
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Figure 15: The helicity amplitude |^l||(s)| 2 in B — > K*£ + £~ at next-to-leading order (solid 
center line) and leading order (dashed). The band for NLO reflects the theoretical uncer- 
tainties from the input parameters. 




Figure 16: The helicity amplitude |-4_i_(s)| 2 in B — > K*£ + £~ at next-to-leading order (solid 
center line) and leading order (dashed). The band for NLO reflects theoretical uncertainties 
from input parameters. 
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2.5 5 7.5 10^2 ^ 2 y- 5 17 - 5 2 " 5 5 7 - 5 1 °^ 2 '^"2'j 15 17-5 



Figure 17: The functions cos <f>u(s) and sin</>n(s) at next-to-leading order (solid center line) 
and leading order (dashed). The band reflects all theoretical uncertainties from parameters 
with most of the uncertainty due to the form factors £i(0). The vertical line at s = 8 GeV 2 
represents the domain of validity of the LEET approach in our case. 




2.5 5 7.5 10.^2^ 2 \ 15 17 - 5 2 " 5 5 7 - 5 10 ^ 2 "^2^ 15 17-5 



Figure 18: The functions cos0^(s) and sin0j_(s) at next-to-leading order (solid center line) 
and leading order (dashed). The band reflects all theoretical uncertainties from parameters 
with most of the uncertainty due to the form factors ^(O).The vertical line at s = 8 GeV 2 
represents the domain of validity of the LEET approach in our case. 
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Figure 19: The phase 0o(s) at next-to-leading order (solid center line) and leading order 
(dashed). 



6 Decay Distributions in B — ► ptvi 

The differential decay rate for B — > p(— > 7r + 7r _ )£z7 £ can be expressed as follows [48-50]: 

G 2 F \V ub \ 2 ^-B( P -> 7T+7T ) 



rfs d cos 9 P dcos9 + d<j) 8(47r) 4 " m| 

{(1- cos ^+) 2 sin 2 ^ p |#+(s)| 2 

+ (1 + cos^+) 2 sin 2 ^ |i/-(s)| 2 
+4 sin 2 # + cos 2 ^ \H (s)\ 2 

—4 sin# + (1 — cos# + ) sin^p cos6 p cos0 H + (s)Ho(s) 

+4 sin# + (l + cos6> + ) sin^ p cos6>p cos H_(s)H (s) 

-2 sin 2 # + sin 2 ^ cos 20 H + (s)H_ . (64) 

The three angles 6+, 9 P and are defined as follows: 9 + is defined by the direction between 
the charged lepton and the recoiling vector meson measured in the W rest frame, the polar 
angle 9 P is defined by the directions of the 7r + (or 7r~) and the vector meson in the parent 
meson's rest frame, and the azimuthal angle is the angle between the two planes, defined 
by the momenta of 7r + 7r~ and the lepton pair £p. 

The helicity amplitudes can in turn be related to the two axial- vector form factors, Ai (s) 
and A 2 (s), and the vector form factor, V(s), which appear in the hadronic current 



H ± (s) = (m B + m p ) A 1 (s) T 2 ^ V(s) , (65) 

tub + m p 



1 



2m p yfs 



A 



H (s) = -= (m 2 B -m 2 -s)(m B + m p ) Ai(s) - 4 ■ A 2 (s) . (66) 



rriB + m. 



Using Eqs.(18),(19) and (20) in Eqs. (66) and (65), we obtain the helicity amplitudes in 
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the large energy Limit: 
H±(s) 
H (s) 



Vx~ 

m B 



Us) , 



m B m p t/s 
2 A 



m B E p ^fs 



m B E p (m| — m 2 — s) — 2X 



(67) 
(68) 



We give below the double differential decay rate (Dalitz distribution) for B — > p(— > ir + n~)lis 
in the variables (s, 0), (s, cos + ) and (s, cos 6> p ), giving also the expressions for the individual 
contributions from the Helicity-0 and Helicity — 1 amplitudes for the Dalitz distribution in 

(s, cos# + ): 



d 2 B 

d<ft ds 



r B 



G 2 F s VX 
192m|7r 4 



\V ub \ 2 (B(p^n + n-)) , 



(69) 



d 2 B 



\H-\ 2 



dcos9 + ds 
d 2 Bu 



= T B 



{\H (s)\ 2 + \H_(s)\ 2 + \H + (s)\ 2 - cos 20 H + (s)} , 

G|S ^|K 6 r(£(p-7r+0){(l + costf + ) 2 \H4s)\ 2 } , 



256m 3 B ir 3 



>|#o| 2 



dcos9 + ds 

d 2 B 
dcos9 + ds 



= r B 



Gl S yf\ 



= r B 



256m| 



G% S V\ 

256mi7T 3 



-|K6| 2 (^(p^vr + 7r-)){2 sin 2 # + \H (s)\ 2 } , 
\V ub \ 2 (B(p^7T + n-)) , 



d 2 B 



dcos9 p ds 



1 2 sin 2 # + \H (s)\ 2 + (l- cos9 + ) 2 \H + (s)\ 2 + (1 + cos# + ) 2 |#-(s)| 2 } 

(70) 

(71) 



d 2 B lHol 2 d 2 B\ H+ \2 d 2 B\ H _\2 



+ 



+ 



dcos9 + ds dcos9 + ds dcos9 + ds ' 
r B G t S ^ \Vu b \ 2 (B(p^^n-)) 



{2 cos 2 9 p \H (s)\ 2 + sin 2 9 p (|iJ + ( S )| 2 + |iJ_( S )| 2 ) } . 



In Figs. (20), (21), (22) and (23), we show, respectively, the Dalitz distributions (d 2 B\ H _\2 /dcos9 + ds), 
(d 2 B\ Ho \ 2 /d cos 9 '+ ds), (d 2 B/dcos9 + ds) and (d 2 B/dcos9 p ds). 

Integrating out the angle 6+, 9 P and from Eq. (64), we obtain the total branching decay 
rate: 



dB G 2 F s y/X 

~ds~ ~ TB 96m|7r 4 



\V ub \ 2 {B{p - n+n-)){\H (s)\ 2 + \H + (s)\ 2 + \H4s)\ 2 } (72) 

dB\H \ 2 dB\ H+ \2 dB\ H _\2 
ds ds ds 

Just as in the decay B — > K*l + l~ , the contribution from the |i/ + (s)| 2 is negligible, and we 
do not show it here. The contributions from the \H_(s)\ 2 , \H (s)\ 2 and the total are shown 
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in Figs. (24), (25) and (26), respectively. The impact of the NLO correction on the various 
branching ratios in B — > p£v^ is less significant than in the B — * K*£ + £~ decay, reflecting 
the absence of the penguin-based amplitudes in the former decay. 

Concerning the B — > p£i>e form factors, one has to consider the SU(3)-breaking effects in 
relating them to the corresponding form factors in B — > K*£ + £~ . For the form factors in full 
QCD, they have been evaluated within the Light-cone QCD sum-rule approach [19,51,52], 
Lattice-QCD [17], and in quark models of the more recent vintage [53,54]. Based on the 
Light-cone QCD sum-rule approach, we estimate the SU(3)-breaking in the ratio of the 
LEET-based form factors at s = as 

Oo) 

C m3) = -$— = 1.3 ±0.06. (73) 

S±,||V U J 

With this and ^ ^(0) from Table 3, we obtain: 

£[ p) (0) = 0.22 ±0.04 . (74) 

The results for the various helicity amplitudes H + (0), H_(0), y/sHo(0), and the ratio 
H_(0)/H + (0) in the decay B — > plv^ are given in Table 5, and compared with some other 
estimates of the same in the literature. We recall that in our approach both ^ ^(0) and 
£^(0) are smaller due to the normalization of the former from data on B — > K*j decay 
rate. This is reflected in the smaller values of the helicity amplitudes in this work compared 
to the other approaches. 

To get the B — > p £ vn form factors at s ^ 0, we use the same parametrization as the 
one for the B — > K*£ + £~ form factors: 



«a« - • (75) 



This parametrization has been used in calculating the differential branching ratios shown 
for the restricted region s < 7 GeV 2 in Figs. (24), (25) and (26). While we do not insist on 
using our approach to cover the small i? p -region (and, hence higher values of s) in the decay 
B — > p£u£, but for the sake of comparison with recent data and some existing estimates 
of the various helicity amplitudes in the literature, we use the s-dependence in (75) to 
estimate H^s) in the entire kinematic region. Numerical results from our work are shown 
in the first row of Table 6 for s = s max - The equality H + (s niax ) = if_(s max ) = H (s max ) 
is imposed by the kinematics of the decay B — > pivg, but we note that our results are 
numerically smaller than the ones following from other form factor models shown in this 
table. This is again to be traced back to our normalization of the function ^ (0) and (75). 
In Table 7, we present various reduced partial widths for B — > p£u^ for the longitudinal 
part Tl = To, transverse part = T_ + r+, the total reduced width r^tai, factoring 
out the CKM matrix element |Kb| 2 , and the ratio T L /T T , and compare them with the 
corresponding estimates in the literature. Our values for the total reduced decay widths 
are smaller than the other five shown in Table 7. For the form factor model in Ref. [53], 
we show a detailed comparison and note that our estimate of (derived from LEET and 
data) is significantly smaller than in this model, but our T L , with input from the Light-cone 
QCD sum rules, is comparable to the one in Ref. [53]. This has the consequence that in 
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our approach Fl and Ft are about equal, as opposed to theirs where F^ is significantly 
smaller than Ft, which will be tested in future experiments. 

Finally, we note that the smaller values of the form factors in the LEET-based approach 
used in this work lead to estimates of \V u b\ from the measured exclusive decay rate for 
B — > pivi> [56], which are in agreement with the PDG 2004 average [57] \V u b\ = (3.67 ± 
0.47) x 10~ 3 . This can also be seen as follows: Using the central values of \V ub \ and the 
lifetime r^o = 1.536 ps from PDG 2004, and the reduced width in the LEET approach 
TTotai = 10.49 ps -1 from Table 7, we get B(B — > ptvi) = 2.18 x 10~ 4 , in comfortable 
agreement with the measured branching ratio by the BABAR collaboration [56] B(B — > 
plvi) = (2.14 ± 0.21 ± 0.48) x 10~ 4 . In contrast, the allowed range of \V u b\ obtained 
recently by the BABAR collaboration [56] using the form factors from the Light-cone QCD 
sum-rules [52] is |K&| = (2.82±0.18±0.30±0.18) x 10~ 3 , with similar values for the quark- 
model based form factors [55]. This is significantly smaller than the PDG 2004 average [57] 
|T4&| = (3.67 ± 0.47) x 10~ 3 . Thus, it is evident that, unless the corrections to the LEET- 
based approach are very significant, data for the radiative decays B — > (K*,p, but 
also for the semileptonic decay B — > pive, favour smaller form factors than obtained in the 
QCD sum rule approaches or the quark models. 





H + (0) 


#-(0) 


H_ (0) /H + (0) 


V~sH Q (0) 


This Work 


0.048 


2.274 


47.002 


7.459 


Ref. [53] 


0.256 


3.253 


12.714 


8.737 


QM [54] 


0.163 


2.942 


18.014 


8.202 


LCSR [19] 


0.055 


3.103 


56.767 


10.158 


Lattice [17] 


0.0898 


4.224 


47.020 


10.914 



Table 5: Various helicity amplitudes for B — > piv^ calculated at s = 0, using different 
models. 





H+ (^max) 


H— (Smax) 


Hq (Smax) 


\/~sHq (s max ) 


This Work 


1.919 


1.919 


1.919 


8.656 


Ref. [53] 


2.362 


2.362 


2.362 


10.652 


QM [54] 


2.952 


2.952 


2.952 


13.313 


LCSR [19] 


2.782 


2.782 


2.782 


12.549 


Lattice [17] 


2.826 


2.826 


2.826 


12.744 



Table 6: Various helicity amplitudes for B — > piv^ calculated at s = s maX ) using different 
models. 
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Figure 24: The dilepton invariant mass distribution dBinw/ds for B — ► piv at next-to- 
leading order (solid center line) and leading order (dashed). The band reflects theoretical 
uncertainties from input parameters. 
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Figure 25: The dilepton invariant mass distribution dB\n \ 2 /ds for B — ► piv at next-to- 
leading order (solid center line) and leading order (dashed). The band reflects theoretical 
uncertainties from input parameters. 
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Figure 26: The dilepton invariant mass distribution for B — > piv at next-to-leading order 
(solid center line) and leading order (dashed). The band reflects theoretical uncertainties 
from input parameters. 





r_ 


r+ 


r L = r 


r T = r_ + r + 


TTotal 


Vl/Tt 


This Work 


4.978 


0.084 


5.423 


5.062 


10.49 


1.07 


Ref. [53] 


9.662 


0.059 


6.095 


9.721 


15.82 


0.63 


QM [54] 










15.8±2.3 


0.88±0.08 


LCSR [19] 










13.5±4.0 


0.52±0.08 


Lattice [17] 










16.533 


0-801^ 


ISGW2 [55] 










14.2 


0.3 



Table 7: Reduced partial and total decay rates for B — > ptvi in units of ps _1 , the ratio 
Ti/rr, obtained in our approach, and using different models for form factors. To get the 
decay widths, one has to multiply the entries in this table with [V^l 2 . 



7 Determination of |\4&|/|Vt s | from B —> plug and 
B -> K*£ + £~ Decays 

The measurement of exclusive B — > piv decays is one of the major goals of B physics. 
It provides a good tool for the extraction of \V u b\, provided the form factors can be either 
measured precisely or calculated from first principles, such as the lattice-QCD framework. 
To reduce the non-perturbative uncertainty in the extraction of V u b, we propose to study 
the ratios of the differential decay rates in B — > plvi and B — > K*£ + £~ involving definite 
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Figure 27: The Ratio R-(s) with three indicated values of the CKM ratio R b = 
VjfeV^*|. The bands reflect the theoretical uncertainty from (su(3) = 1-3 ± 0.06 and 
^* } (0) = 0.28 ±0.04. 



helicity states. These s-dependent ratios Ri(s), {i = 0, —1, +1) are defined as follows: 

, dTfr K * l+l ~/ds , 
R^s) = _^ I— 76 

The ratio R-{s) suggests itself as the most interesting one, as the form factor dependence 
essentially cancels (in the SU(3)-symmetry limit). From this, one can measure the ratio 
| | / 1 V^tb | . In Fig. (27), we plot R-(s) for three representative values of the CKM ratio 
Rb = \V ub \/\V tb V t *\ = \V ub \/\V cb \ = 0.08, 0.094, and 0.11. We also show the ratio R (s), 
where the form factor dependence does not cancel. For the LEET form factors used here, the 
compounded theoretical uncertainty is shown by the shaded regions. This figure suggests 
that high statistics experiments may be able to determine the CKM-ratio from measuring 
R-{s) (and Rq(s)) at a competitive level compared to the other methods en vogue in 
experimental studies. 
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Figure 28: The Ratio Rq(s) with three indicated values of the CKM ratio i?& = | V^&|/| Vti,V^*|. 
The bands reflect the theoretical uncertainty from (su{3) = 1-3 ± 0.06 and £± \o) = 
0.28 ± 0.04. 
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8 The Ratios i?_(s) and Rq(s) as Probes of New Physics 

in B -> K* l+r 



In order to look for new physics in B — > K*£ + £~ , we propose to study the ratios Ro(s) 
and R-(s), introduced in the previous section. As well known, new physics can distort the 
dilepton invariant mass spectrum and the forward-backward asymmetry in a non-trivial 
way. 

To illustrate generic SUSY effects in B — > K*£ + £~, we note that the Wilson coefficients 
Cf ff , Cf ff , Co, and C w receive additional contributions from the supersymmetric particles. 
We incorporate these effects by assuming that the ratios of the Wilson coefficients in these 
theories and the SM deviate from 1. These ratios for k — 7, 8, 9, 10 are defined as follows: 



They depend on the renormalization scale (except for C w ), for which we take /x = m bjPole . 
For the sake of illustration, we use representative values for the large-tan (3 SUGRA model, 
in which the ratios r-j and r§ actually change their signs. The supersymmetric effects on 
the other two Wilson coefficients Cg and C\q are generally small in the SUGRA models, 
leaving r 9 and no practically unchanged from their SM value. To be specific, we take * 



In Figs. (29) and (30), we present a comparative study of the SM and SUGRA partial 
distribution for H and H Q , respectively. In doing this, we also show the attendant theo- 
retical uncertainties for the SM, worked out in the LEET approach in this paper. For these 
distributions, we have used the form factors from [18] with the SU(3)-symmetry breaking 
parameter taken in the range Csu(3) — 1-3 ± 0.06. Figs. (29) and (30) illustrate clearly 
that despite non-perturbative uncertainties, it is possible, in principle, in the low s region 
to distinguish between the SM and a SUGRA-type models, provided the ratios r k differ 
sufficiently from 1. 

9 Summary and Concluding Remarks 

Summarizing briefly our results, we have reported an 0(a s )-improved analysis of the 
various helicity components in the decays B — > K*£ + £~ and B — > p£u e , carried out in 
the context of the Large-Energy-Effective-Theory. The underlying symmetries in the large 
energy limit lead to an enormous simplification as they reduce the number of independent 
form factors in these decays. The LEET-symmetries are broken by QCD corrections, and we 
have calculated the helicity components implementing the 0(a s ) corrections. The results 
presented here make use of the form factors calculated in the QCD sum rule approach. The 
LEET form factor £±*(0) is constrained by current data on B — > K*^y. As the theoretical 
analysis is restricted to the lower part of the dilepton invariant mass region in B — > K*£ + £~ , 
typically s < 8 GeV 2 , errors in this form factor are not expected to severely limit theoretical 

*We thank Enrico Lunghi for providing us with these numbers. 




(77) 



= -1.2, r 8 = -1, r 9 = 1.03, r 10 = 1.0 . 



(78) 
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Figure 29: The Ratio R_(s) with \V ub \/\V tb V* s \ = 0.094 in the Standard Model and in 
SUGRA, with (r 7 , r 8 ) = (—1.2, —1), (su(3) — 1-3 an d ^± \o) = 0.28 represented, respec- 
tively, by the shaded area and the solid curve. The shaded area depicts the theoretical 
uncertainty ( SU (3) = 1.3 ± 0.06 and ^* } (0) = 0.28 ± 0.04. 



precision. This implies that distributions involving the H^(s) helicity component can be 
calculated reliably. Precise measurements of the two LEET form factors £±( s ) anc ^ £\\( s ) 
in the decays B — > ptvi can be used to largely reduce the residual model dependence. 
With the assumed form factors, we have worked out a number of single and double (Dalitz) 
distributions in B —>■ plug, which need to be confronted with data. An analysis of the decays 
B — > K*£ + £~ is also carried out in the so-called transversity basis. We have compared the 
LEET-based amplitudes in this basis with the data currently available on B —>■ K*J/i/)(—> 
£ + £~) and find that the short- distance based transversity amplitudes are very similar to 
their long-distance counterparts. We also show the 0(a s ) effects on the forward-backward 
asymmetry, confirming essentially the earlier work of Beneke, Feldmann and Seidel [13]. 
Combining the analysis of the decay modes B — > K*£ + £~ and B — > p£vi, we show that 
the ratios of differential decay rates involving definite helicity states, R-(s) and Ro(s), 
can be used for testing the SM precisely. We work out the dependence of these ratios 
on the CKM matrix elements |T4b|/|^ts|- We have also analyzed possible effects on these 
ratios from New Physics contributions, exemplified by representative values for the effective 
Wilson coefficients in the large-tan (3 SUGRA models. The main thrust of this paper lies, 
however, on showing that the currently prevailing theoretical uncertainties on the SM 
distributions in B — > K*£ + £~ can be largely reduced by using the LEET approach and 
data on5-^ K #r y and B — > p£i>£ decays. Finally, we remark that the current experimental 
limits on B -> (X s ,K*)£ + £~ decays (and the observed B -> K£+£r decay) [39-42] are 
already probing the SM-sensitivity. With the integrated luminosities over the next couple 
of years at the B factories, the helicity analysis in B — > p£v^ and B — > K*£ + £~ decays 
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Figure 30: The Ratio Rq{s)/ with \V ub \/\V tb V* s \ = 0.094 in the Standard Model and in 
SUGRA, with (r 7 , r 8 ) = (—1.2, —1), (su(3) — 1-3 an d ^(0) = 0.28 represented, respec- 
tively, by the shaded area and the solid curve. The shaded area depicts the theoretical 
uncertainty ( SU (3) = 1.3 ± 0.06 and ^* } (0) = 0.28 ± 0.04. 



presented here can be carried out experimentally. 
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